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THE DISPLACEMENTS IN AN ELASTIC HALF-SPACE DUE 
TO A MOVING CONCENTRATED NORMAL LOAD 

By Donald L. Lansing 
Langley Research Center 

SUMMARY 

Equations are derived for the steady- state displacements produced at any point 
within an elastic half -space by a concentrated normal load in uniform motion over the 
surface. The results are valid for all depths, load speeds, and elastic constants. The 
equations for the displacements produced on the surface and in the vertical plane con- 
taining the line of motion of the load for the special case of Poisson's ratio equal to 1/4 
are obtained from the general results and are expressed in terms of elementary func- 
tions and elliptic integrals. Some numerical calculations based upon these equations are 
presented and discussed for a wide range of load speeds. In addition, the complete time- 
dependent equations for the horizontal displacements produced at the surface, including 
the transients excited by the initial application of the load, are derived. 

INTRODUCTION 

Supersonic aircraft, air bomb bursts, and surface vehicles produce forces which 
act upon and move over the earth's surface. As a result of the earth's elasticity, these 
loads are transmitted through the earth's outer layers and over the ground surface and 
thereby produce disturbances at some distance from the area directly in contact with the 
load. These disturbances can set off a buried mine, can shake building occupants, and, if 
sufficiently severe, can produce damage to structures. The practical need for estimating 
the magnitudes of the ground motions and stresses induced by moving surface loads as 
well as the interesting theoretical aspects of the problem has produced over the past few 
years a number of papers on the fundamental dynamic phenomena involved. (See refs. 1 
to 5.) 

In order to investigate the displacements and stresses produced in the earth by a 
moving surface load by using the classical linear theory of elasticity, it is necessary to 
idealize the actual situation to some extent. A convenient starting point to the problem 
of calculating these disturbances is to assume that they are produced on a uniform semi- 
infinite elastic half space by a concentrated normal surface load moving in a straight line 
at constant speed. No generality is lost in assuming a concentrated load since more 



general load distributions may then be built up by familiar superposition principles. The 
assumptions of a uniform medium with a plane surface and rectilinear motion limit the 
usefulness of the results but these conditions are fulfilled, or nearly so, in some practi- 
cal problems and are sufficient for gaining an insight into the basic phenomena to be 
expected in more realistic situations. 

A theoretical analysis of even this simplified problem is made difficult, and chal- 
lenging, by the fact that there are three fundamental wave speeds in an elastic half -space. 
In decreasing order of magnitude, these speeds are the speeds of propagation of longitu- 
dinal waves, shear waves, and surface waves. The solution takes different forms and 
the details of an analysis vary depending upon the speed of the load relative to these three 
wave speeds. In this paper the load speed is referred to as "subseismic,” "super- 
seismic," and "transeismic" when it is less than the shear wave speed, greater than the 
longitudinal wave speed, and intermediate between these two speeds. Furthermore, one 
may decide to treat the complete time-dependent problem and account for the starting 
transients excited by the initial application of the load or, in the interests of further sim- 
plicity, one may treat just the steady-state case and obtain a description valid only after 
the load has been moving for some time. These factors and other considerations lead to 
a wide variety of partial problems which are gradually being solved piecemeal as the 
need for further information and theoretical interest in the basic problem develops. 

The two-dimensional moving-line load was treated in the steady-state case in ref- 
erence 1 at subseismic speeds and in reference 2 where the displacements and stresses 
were derived for all load speeds. These results were extended at subseismic speeds in 
reference 3 to include in the stresses the transient effects excited by the initial applica- 
tion of the line load. The steady-state displacements produced by a moving-point load at 
low subseismic speeds were derived in reference 4. The transient solution for the sur- 
face displacements was treated at all load speeds in reference 5. No work appears to be 
available on the subsurface stresses produced by the moving-point load. 

The present paper is a continuation of the study of the three-dimensional case. 

Two aspects of this problem are considered. First, equations are derived for the 
steady-state displacements produced within an elastic half-space by a concentrated nor- 
mal load moving at constant speed in a straight line over the surface. The results, 
obtained by Fourier transform methods, are valid for any depth, load speed, or elastic 
constants and could be used for an evaluation of the subsurface stresses and strains. 

The general equations for the displacements, which are in the form of double integrals 
with infinite limits, are then reduced to single integrals with finite limits. It is further 
shown that in the special case for Poisson’s ratio equal to 1/4, the surface displace- 
ments can be expressed in terms of elementary functions and elliptic integrals. Some 
numerical calculations of the steady- state displacements at the surface and in the 
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vertical plane containing the line of motion of the load are presented and discussed. Sec- 
ond, equations are derived for the horizontal surface displacements which include the 
transients excited by the initial application of the load. These results, obtained by super- 
imposing the response of the surface to an impulse, can also be expressed in terms of 
elliptic integrals and contain the steady-state solution at the surface as a special limiting 
case. 


SYMBOLS 


a,b,c 

A(0),B(e) 

c(e),D(e) 

B I =\|l -M z 2 

Bt = \Jl - M t 2 



H 

C R 

c t 


variables 

quantities defined in equation (46) 


velocity of longitudinal compression waves, 


velocity of Rayleigh waves 
velocity of transverse shear waves, 



+ 2/i 1 

P o 


C+(a,b,c) 

C“(a,b,c) 

f,g 


a^ + (b - c)2 + (b + c)2 

a a 

a^ + (b - c)2 a 2 + (b + c)2 

functions 


H(x) 


unit step function, +1 when x > 0 and 0 when x < 0 



k 

K(k) 


m r = I- 

R C R 
v 

M t = £ 


modulus of elliptic integrals 


complete elliptic integral of the first kind 



de 


- k 2 sin 2 0 


n 1 ( 1 //),n 2 (^) 

P 



r,^ 

r,<Pj 


quantities defined in equations (61) and (62) 
magnitude of concentrated load 


polar coordinates 

radial displacement produced on surface by an impulse 


R(t) = ^R 5 


Sj = \j\2 + ju 2 - X 2 M^ 2 - iVeX 
^ = ]p + ju 2 - X 2 M t 2 - iVeX 


S + (a,b,c) = b - c ■ . + b + c 


a z + (b - c) 2 a 2 + (b + c) 


sgn(x) 


sign function, 

M 
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t 


time 


T((p) = ^1 - M^sin 2 ^ ^1 
T(x,y) = (r 2 - I M t 2 y 2 ) (r 2 


Q!j 2 M^ 2 sin 2 (^ ^1 - a^M 2 sin^ cpj 


Q!l 2 M t 2 y 2 j ( r 2 - Q<2 2 Mt 2 y 2 


u,v,w 


U’jV'jW' 


U M jV M ,W M 


u l> u 2> u 3 

vi,v 2 ,v 3 

w l> w 2> w 3 


Cartesian components of displacements 

contribution of Rayleigh poles to displacements (see eqs. (27) 
and (49)) 

integrals in displacements (see eqs. (27) to (47)) 
quantities defined in equations (37) to (45) 


UpU y 


polar components of steady-state surface displacements 


u r ,u e 

vv 


v 


x,y,z 

x = x - Vt 


“2 


2 = 


3 - \f3 
4 


polar components of transient surface displacements 

polar components of steady-state displacements in vertical 
plane y = 0 

velocity of moving load 

Cartesian coordinates along X-, Y-, and Z-axes, respectively 
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h = 

y l> y 2’ y 3 constants 

8(x) Dirac delta function 

r i 2 

A(A,m) = 2(a 2 + ii 2 ) - M t 2 A 2 - 2ieVAj - 4S^St(x 2 + ii 2 - ieVx) = A(A,m,VA) 

A(A,fi,a>) quantity defined in equation (21) 

e a small positive constant 

2M t 2 sin 2 <p 

r? = 2 2 

3 - M t z sin 

= sin"-*- 




0t = sin ''(^) 


K = k 


\,ll,w 

transform variables 

A',m’ 

Lame constants 
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MA 2 


Rayleigh poles 


n(o! } k) 

P 0 

p,v 

°xx>°yy >°zz 



T xy> T xz> T yz 

?>X 


branch points 

complete elliptic integral of the third kind, 



density of medium composing the half space 
polar coordinates 
normal stresses 

shear stresses 
stress functions 
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denotes an integral taken in sense of a Cauchy principal value 




three-dimensional Laplacian operator 
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STATEMENT OF THE PROBLEM 


Let a concentrated load of magnitude P be applied normal to the surface of a 
semi-infinite elastic half -space and be moved along a straight line at a constant velocity 
V as shown in figure 1. Introduce a right-handed system of rectangular Cartesian coor- 
dinate axes X, Y, and Z fixed in the 
medium with the positive X-axis coinci- 
dent with the line of motion and oriented 
in the direction of motion, and with the 
positive Z-axis pointing into the half 
space. The sign of the load P will be 
taken as positive when it acts in the 
Figure 1.- Coordinate system and components of displacement. direction of the positive Z-axis, that is, 

when the load places the adjacent surface 

in compression. Let u, v, and w be the displacements of any point of the medium 
parallel to the X-, Y-, and Z-axes, respectively. The problem to be solved in this paper 
is the determination of the displacements for any value of the load velocity V. 

Two aspects of this problem are considered. First, it is assumed that the load 
has been moving for a considerable length of time so that all transient disturbances 
produced at the initial point of application have died out and have left only a steady pat- 
tern of displacements with respect to an observer moving with the load. By using 
Fourier transform techniques or, equivalently, superpositions of plane waves, expres- 
sions are obtained from which the displacements at any depth may be calculated. Second, 
the complete time-dependent solution for the horizontal displacements produced at the 
surface of the medium is discussed. In this case, the displacements are obtained by a 
suitable integration in time of the response of the surface to an impulse and include the 
transients excited by the initial application of the load. 

STEADY-STATE DISPLACEMENTS 

The determination of the steady-state displacements requires finding a solution of 
the equations of motion (ref. 6, ch. 5) 

p T V 2 u + (V + M’)V(V . u) = p Q (1) 

dr 

which satisfies at the surface z = 0 the conditions that the normal stress must equal 
-P6(y) S(x - Vt) and the shear stresses must vanish. In these equations u is the dis- 
placement vector ui+vj+wk, p 0 is the mass density of the medium (assumed to be 



8 



constant), X' and p' are the Lame constants of the medium, and S(y) and 5(x - Vt) 
are Dirac delta functions. A solution to this problem is obtained by expressing the dis- 
placements and stresses in terms of two stress functions and then representing the 
stress functions as appropriate superpositions of plane waves to satisfy the surface 
conditions. 


As shown 


Stress Functions and the 
by Lamb (ref. 7), the equations 


= + 

9 2 X 

8x 

ax az 

= + 

9 2 X 

9y 

9y 9z 

= *<£.+ 


9Z 

8z 2 


Radiation Condition 

of motion (eq. (1)) are satisfied by letting 


JL£!z 

c t 2 at 2 
1 J 


( 2 ) 


provided that the two stress functions cp and y are solutions of the three-dimensional 
wave equations 


V 2 cp = 


1 9 2 <p 
c I 2 dt 2 


(3) 


and 


v 2 x 


l a 2 x 
c t 2 at 2 


(4) 


in which C ; 2 = -i-(X + 2ju') and c+ 2 = — p*. In terms of these stress functions the nor- 
Po Po 

mal stresses become 


cp at 2 lax 2 az ax 2 ) 


) 2 7n a 3 v 


o YY =*L d h + 

Ci 2 8t 2 \8y 2 8z 9y 2 J 
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CT ZZ - 


c, 2 8t 2 


9 ,.9 2 ^ ® 3 X 

+ 2/i — £ + — 

\8z 2 9z‘ 


2 Po 


s3 X r 
8z 8t 5 


(5) 


and the shear stresses become 


T xy - T yx 


= ^2L_) 

^ \8x 8y 8x 8y 8z J 


T xz = T zx = 2fi'l-?-JL + d * 


T yz = T zy = 2 m’ 



( 6 ) 

( 7 ) 


where is the speed of propagation of longitudinal compression waves which are 
referred to as P-waves and c^ is the speed of propagation of transverse shear waves 
referred to as S-waves. (See ref. 6, ch. 13.) Both of these waves are three-dimensional 
body waves which spread through the interior of the medium and decay as the reciprocal 
of the distance from the source. 

In addition to these two types of waves, the presence of a free surface on the solid 
gives rise to a third type of wave known as a Rayleigh wave or R-wave. (See ref. 6, 
ch. 13.) This wave is a surface wave similar to a gravity wave on water. It involves 
both compression and shear. The Rayleigh wave is essentially two dimensional; its 
strength falls off rapidly with depth and the energy it carries is confined to a thin layer 
near the surface. 

As will be seen subsequently, the three wave speeds cr, c^, and c^ play an 
important role in determining the nature of the displacement pattern. The P-wave is the 
fastest traveling wave. The Rayleigh wave travels with a speed Cjj which is slightly 
less than the shear wave speed. 

Besides satisfying the imposed boundary conditions, the solution of the present 
problem must also satisfy the radiation condition which is the physical requirement that 
only outgoing waves are observed at a great distance from the applied load. (See ref. 8.) 
As discussed in reference 9 this condition may be taken into account by adding the terms 

e ^ and e ^ (with e positive) to the right-hand side of equations (3) and (4), 

respectively, and considering the original problem as the limiting case e approaching 

0. 
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Accordingly, in place of equations (3) and (4), it is advantageous to deal with the 
equations 


V 2 <p = 


1 9 2 ^ 

c l 2 at2 


+ 6 



and 


V 2 X = 


J-^Z +e 

c t 2 at 2 


9t 


( 8 ) 


(9) 


These equations possess plane-wave solutions of the form 


and 


where 


- = giAx+ijuy-z^-icut 
— _ iAx+i/iy-zffy-icot 

A ^ 


^=|x 2 + M 2.^|. icoe 

and 

= jx 2 + m 2 - - i^e 

Since equations (8) and (9) are linear, these simple solutions may be superimposed by 
summation or integration to construct more general solutions. 

For present purposes it is convenient to assume the following forms for the stress 
functions: 


J +°° n +°o p 

-OO - oo J - 


<p 


f(A,/i,co)e 


iAx+i jiy - icot - z Sfy 


dA dj u. du> 


( 10 ) 


p +°o p +oo p +oo 

d — OO d_oo d_ 0O 


. . iAx+ijuy-icot-z^. , 

g(A,/i,u>)e dA d/i do> 


(ID 


In equations (10) and (11), f(A,p.,a>) and g(A,/i,u>) are arbitrary functions. Once these 
functions have been determined so that the boundary conditions at the surface are satis- 
fied, the displacements may be found from equation (2). The necessary steps in finding 
f and g and then u, v, and w are now carried out. 
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Formal Solution for Displacements 
The boundary conditions to be satisfied at z = 0 are 

or zz = -P5(y)6(x - Vt) 


and 


(12) 


T xz - T yz - 0 


(13) 


From equations (6) and (7) it is seen that equation (13) can be satisfied by setting 


+ 1 £Xo 

ct 2 ae 


} 2 — 


(az ' 8z 2 ) n ,2 a+2 


(14) 


at z = 0. Equations (14), (5), and (12) give two conditions for determining the unknown 
functions f and g. If the expressions for 7p and y (eqs. (10) and (11) are substi- 
tuted into equations (5) and (14) and z is set equal to zero, one obtains, with the aid of 
equation (12), 


r> +°° p +°° n +°° j 

■P6(y)6(x - Vt) = \ \ \ f ( 2 p 

'- 1 — CO 'J — OO ^—00 \ 


•a 


l 


2 

c l 2 J 


- 2fi t M , gfot 2 + ^2 


e iAx+ipy-iwt dx djj da) 


(15) 


and 


» +°° r* +°° r* + 00 


°= y j j 

— 00 ^ -oo ^ -c 


2fyf - g[2S2t 2 


e iAx+ip.y-icot dA djU ^ 


(16) 


Equations (16) and (15) may be solved for the terms in brackets in the integrands by 
means of the Fourier inversion formulas (ref. 10): 


2%f - g\2Sl t 2 = 0 


(17) 


12 



and 


f 



x '$)- 2nt " , ( n ‘ 2 + $) B 


/-» "f*°0 /-fc +00 -f-oo 

= - \ \ \ 6(y)6(x - Vt)e 

x O J _oo *J-oo -co 

(2tt) 

= - 6(o> - VX) 

(2tt) 2 


iXx-i M y+ia>t dxdydt 


When solved for f and g, equations (17) and (18) give 

2«t 2 +^ 

f = ^^7 ^-6(c-VX) 

(277 r A U>M,U>) 


(18) 


(19) 


g = 


-P' 


2f2, 


(2tt) 2 a ( x >M,w) 


6(u> - VX) 


where P’ = — r and 

M t 


A(X,M,w) = 


2(X 2 + M 2 ) - ^ - 2ieco 
ct 2 


4fyS2t(x2 + jLt 2 - ieoj) 


( 20 ) 


( 21 ) 


The stress functions <p and y are now determined by substituting equations (19) 
and (20) into equations (10) and (11). The integration with respect to co is easily per- 
formed and yields 


and 


<P = ~^2 

477^ 


is : 


+°o /-» +0O 
—oo 


2 St 2 + M t 2 X 2 e iXx+ifxy-zS i dx 

A(X,m) 


(22) 


where 


i.^rr j5l. e^y-^ <,x dM 

477 2 ^-oo A (X,p.) 


x = x - Vt 



(23) 
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and, from equation (21), 


Si 2 = X 2 + IX 2 - M t 2 X 2 - iVeX 
St 2 = X 2 + \L 2 - M t 2 X 2 - iVeX 


A (X,p) = A(X,p,XV) = ^(X 2 + p 2 ) - M t 2 X 2 - 2ieVx] - 4S^(X 2 + is 2 - ieVX) 

The displacements obtained by substituting equations (22) and (23) into equation (2) are 

u<s ’ y ’ z > - B 07 "S^h 2 + M ‘ 2i2)e ' zs ' - 2s ' d " <24 > 

v(x,y,z) = =¥± f + f + ^ Lst 2 + M t 2 X 2 )e“ zS i - 2S Z Ste _zS t]dX d p (25) 

4 2 J-oo^.oo A(x,/i) L J 


w(x,y,z) = f f ^ — r^|(2St 2 + Mt 2 X 2 )e _zS ^ - 2(St 2 + Mt 2 X 2 )e _zS t]dX dp (26) 
4 7T ^ — oo /ij L J 


Since these expressions depend upon t only through the variable x, the displacement 
pattern is steady in time for an observer moving with the load. 

Equations (24) to (26) are the general expressions for the displacements valid for 
any depth, load speed, or elastic constants. In particular, when V = 0, these equations 
reproduce the classical solution of the problem of Boussinesq. (See ref. 6, ch. 8.) In 
addition, it is now a straightforward matter to obtain corresponding integrals for the 
stresses and strains produced within the medium. This phase of the problem will not be 
pursued here and attention will be focused on the displacements only. 

The integrands of equations (24) to (26) are, in general, multiple -valued functions 
because of the presence of the radicals Sj and St and become infinite at the so-called 
Rayleigh poles which occur at a certain pair of simple zeros of A (A,/i). Therefore, the 
value of the radicals must be defined and the manner of integrating around the poles must 
be determined before the integrals can be evaluated numerically or simplified by carrying 
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out some of the integrations in closed form. The proper treatment of these singularities 
which is necessary to obtain a physically acceptable result requires a detailed study of 
the poles and branch points of the integrands and depends upon the values of the elastic 
constants A' and /j' (or equivalently, c^ and c t ) and load velocity. It is in this 
investigation that the positive parameter e plays a significant role. 

In general, it is found that the expressions for the displacements can be written in 
the form 


(u,v,w) = (u',v r ,w') + (u",v”,w") (27) 

The terms u', v', and w' are associated with the Rayleigh poles and must be taken 
into account whenever the load speed exceeds the Rayleigh wave speed. To obtain explicit 
expressions for these quantities, one must determine the zeros of A(A,ju). This deter- 
mination of zeros requires some assumption for the value of the ratio of elastic constants 
A r and p' or for the value of the ratio of the wave speeds c i and q-. The terms 
u", v", and w" are integrals of the same form as those in equations (24), (25), and (26) 
but are taken in the sense of Cauchy principal values. These double integrals with infi- 
nite limits can be expressed in terms of single integrals with finite limits by transforming 
from the "rectangular" coordinate system A ,p to a "polar" coordinate system r,0 and 
then carrying out the r-integration. The results do not depend upon any assumptions for 
the values of the elastic constants or wave speeds. The integrands and limits of the 
remaining 0 integrals take various forms depending on whether V < q, the subseismic 
case; q < V < q, transeismic case; or q < V, the superseismic case. Explicit 
expressions for the quantities u", v", and w" in the three different cases are given 
in the next section. 


Single Finite Integrals for u", v", and w" 

In this section single integrals with finite limits are presented for the quantities 
u", v", and w". The detailed steps leading from equations (24), (25), and (26) to these 

results are given in the appendix. In the following equations 6^ = sin - ^ and 

indicates that an integral is to be understood in the 

The parameter e does not appear in the final equa- 
tions; its only purpose is to determine the manner of handling the branch points and 
poles. Once this is done, as in the appendix, e is set equal to zero. 


0 7 = sin"! — The symbol X d0 
sense of a Cauchy principal value. 
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The desired expressions for single integrals with finite limits are: 
Subseismic case; 0 < V < c^, < Mt < 1: 

p» p 77 /^ 

u"(x,y,z) = — x ui(x,y,z;0)d0 

7 7 2 *'0 

p» 7 t/2 

v"(x,y,z) = -^ X v 1 (x,y,z;0)d0 

7 T 2 J 0 

pt p 77 /*^ 

w"(x,y,z) = — X w 1 (x,y,z;0)d0 

7T 2 *''0 


Transeismic case; c^- < V < c^, < 1 < M^: 


pt p$ t p 77 / 2 _ 

i"(x,y,z) = i-r X u 1 (x,y,z;0)d0 + X u 2 (x,y,z;0)d0 

77 2 ‘'O 

- 


__ pt p &t _ p 77 /^ _ 

v"(x,y,z) = ^ X v 1 (x,y,z;6)d0 + ¥ v 2 (x,y,z;S)dS 

7T 2 

pt p 9-^ _ p 77 /2 __ 

/”(x,y,z) = -=■ X wi(x,y,z;0)d0 + X w 2 (x,y,z;0)d0 
7r 2 J 0 J St- 


Super seismic case; cj; < c^ < V, 1 < < Mj: 


p, />0t p07 _ p 77 / 2 _ 

u"(x,y,z) = X u 1 (x,y,z;0)d0 + X u 2 (x,y,z;0)dS + X u 3 (x,y,z;S)dS (34) 

7T 2 *^0 J 0t J 9l 

x>, P (9f _ p 0£ _ p7r/2 _ 

v"(x,y,z) = X v 1 (x,y,z;0)d0 + X v 2 (x,y,z;0)dS + * v 3 (x,y,z;0)dS (35) 

77 2 *'0 J 9t J 9i 

p t pSt _ C®1 — p 77 / 2 _ 

w"(x,y,z) = ^ x w 1 (x,y,z;S)dS + X w 2 (x,y,z;0)d0 + X w 3 (x,y,z;0)dS (36) 

7T 2 J o J e t J 9i 
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in which 


ui(x,y,z;0) = - sin S + (z\fl - M z 2 sin 2 0,x sin 0,y cos 0 

A 2 (0) - 4B(0) 2 \ « 


- B(0)S + (z\ 1 - Mt 2 sin 2 0,x sin 0, y cos 0 


vi(S ’ y ’ z;e) = afiiiif f s+ ^ ] - M ^ slnZe ’ y _ 

- B(9)S + fz\Jl - Mt 2 sin 2 0,y cos 0,x sin 0 


Mp-S+Wi - M z 2 sin 2 0,y cos 0,x sin 0 


\/l - M? 2 sin 2 0 a (a) ( I 2 o - 

t (x,y,z;0) = -*-5 C+ zUl - sin 2 0,x sin 0,y cos 0 

A 2 (0) - 4B(0) 2 \ v > 


C + (z\ll - M t 2 sin 2 0,x sin 0,y cos 0 


sin 0 


A 4 (0) + 16C 2 (0) 2 


A 3 (0) 


S + (z\l - M^ 2 sin 2 0,x sin 0,y cos 0 


- 2A(0)C(0)C + ^z\jl - 2 sin 2 0,x sin 0,y cos 0j 
+ 7 r A 2 ( 0) C ( 0) 5^x sin 0 + y cos 0 + z^Mt 2 sin 2 0 - l) 
+ 6^x sin 0 - y cos 0 + z^Mt 2 sin 2 0 - ljj 


8C 2 (0) 


x sin 0 + z\|Mt 2 sin 2 0 - 1 


(x sin 0 + z\|M t 2 sin 2 0 - 1 ) - y 2 cos 2 0 


cos 0 


v 2 (x 5 y,z;0) = 


I A 3 (g) 
A4(0) + 16C 2 (0) ) 2 


cos 0,x sin 0^ 


+ 2A(0)C(0)C 


2 sin 2 0,y cos 0,x sin 0 


- ttA 2 (0)C(0) 


-(-IF 

5^x sin 0 - y cos 0 + z^Mt 2 sin 2 0 - 1 


6 ^x sin 0 + y cos 0 + z^M^Sin 2 0 - 1 


8C 2 (0) 


y cos 0 


x sin 0 + zyM^sin 2 0 - 1^ - y 2 cos 2 0 


w 2 (x,y 


>y> z >p; . a / .V . __o/ -v \ o 


A 4 (0) + 16C 2 (0) 


C + (zul - M£ 2 sin 2 0 3 x sin 0,y cos 0) 


v. 


+ 2A(0)C(0)S + lz\|l - M^ 2 sin 2 0,x sin 0, y cos 0 


7lA 2 (0) 


\ 

6^x sin 0 + y cos 0 + z^Mj- 2 sin 2 0 - 1^ 


+ S^x sin 0 - y cos 0 + z^Mt 2 sin 2 0 - 1 J 
- 8C(0) 


x sin 0 + z^Mt^sin 2 0 - 1 


^x sin 0 + z^Mt 2 sin 2 0 - l) - y 2 cos 2 0 


u 3 (x,y, z ;0) 


v 3 (x,y,z;0) 


w 3 (x,y,z;0) 


sin 0 


A 2 (0) + 4D(0) 


Me) 


x sin 9 + zyM£ 2 sin 2 0 - 1 


^x sin 9 + z(M^) - y 2 cos 2 0 


+ 2D(0) 


x sin 9 + zyMt 2 sin 2 0 - 1 


^x sin 6 + z^M|- 2 sin 2 0 - lj - y 2 cos 2 0 


( 43 ) 


-cos e 


A 2 (0) + 4D(0) 


A(0) ^L^osi 

(x sin 9 + z^Mj 2 sin 2 0 - 1^ - y 2 cos 2 0 


+ 2D(0) 


y cos 0 


(x sin 0 + z^Mt 2 sin 2 0 - 1 

\ 2 

) - y 2 cos 2 0 


(44) 




2 sin 2 0 - 1 


A 2 (0) + 4D(0) 


A(0) 


x sin 0 + z\Mjj 2 sin 2 0 -'1 


^x sin 0 + z^M^siii^9~-~i > j - 


y 2 cos 2 0 


- 2 


x sin 0 + z^Mt 2 sin 2 0 - 1 


^x sin 0 + z^M t 2 sin 2 0 - 1 'j - y 2 cos 2 0 


(45) 
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A(0) = 2 - M t 2 sin2e 


■\ 


and 


B(0) = 

^1 - M t 2 sin 2 0 1 

||l - Mj 2 sin 2 0 

C(0) = 

^Mj- 2 sin 2 0 - 1^ 

jl - M z 2 sin 2 0 

D(0) = 

\|M t 2 sin 2 0 - 1 \ 

Im^ 2 sin 2 0 - 1 

F 


V I 

J 


( 46 ) 


S + (a,b,c) = , b - - c „ + b + c 


a 2 + (b - c) 2 a 2 + (b + c) 


C+(a,b,c) = — s + fi -V 

a 2 + (b - c) 2 a 2 + (b + c) 2 } 


C”(a,b,c) = 


a 2 + (b - c) 2 a 2 + (b + e) 2 


(47) 


When Cpj < V, the integrals of uj, v^, and wi are singular since A 2 (0) - 4B(0) 

C R 

vanishes when 0 = 0pj = sin-1 — which is smaller than either of the upper limits 

or it/2. The integrals of the delta functions in U£, V 2 , and W 2 can be carried out 
analytically by means of the following formula (ref. 11): 


J b f(0)5(g(0))d0 = £ 


f(gj) 

g’(0i)| 


(48) 


in which the summation extends over all the zeros of g (0), the 0j values, which lie 
between the upper and lower limits of integration, a and b. Equations equivalent to 
equations (28) to (30) and (37) to (39) for the subseismic case are given in reference 4. 
However, the terms which must be added to these integrals to obtain the displacements 
when cjj < V < c t are not indicated. 
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Expressions for u', v', and w' When Poisson’s Ratio Is 1/4 

In the remainder of this paper it is assumed that Poisson's ratio is 1/4; thus, 

X' = n’ = 2/5 Young's modulus. The relation between the P-wave speed Cjr and the 
S-wave speed is then thus, the compression waves travel with a speed 

about 70 percent greater than that of the shear waves. The relation between the 

2 

Rayleigh wave speed c R and the S-wave speed is c R = . c t * 0.92^. 

\j3 + y3 

This assumption is reasonably close to the truth for many materials. Moreover, 
the steps involved in obtaining explicit expressions for u’, v’, and w’ in this special 
case indicate the procedure to be followed when other values of X* and ji' are of 
interest. The details are worked out in the appendix; only the final results are given 
here. When the equations are presented, it is convenient to use the abbreviations: 



0R = i/m R 2 - 1 

o 3 + \/3 

“1 2 = — 4 


2 4 


Again it should be noticed that the parameter e does not appear in the final results. 

The additional terms arising from the Rayleigh poles when Poisson's ratio is 1/4 

are: 


u ’(x,y,z) 


v’(x,y,z) 


w'(x,y,z) 
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in which C+(a,b,c), C~(a,b,c), and S + (a,b,c) are the combinations of elementary func- 

tions defined in equation (47). These terms must be included in equation (27) for the 
displacements when Mr > 1 and are to be ignored when Mr <1. In particular, at the 
surface (z = 0) 

u'(x,y,0) = v’(x,y,0) = 0 

and 


w'(x,y,0) 


P'M t \/5 + 3# x 
16 tt/3 r x 2 - /3 R 2 y 2 


(50) 


Equation (49) and the expressions for the stresses associated with these "displace- 
ments" can be used without the additional and more complicated terms u", v", and w" 
to estimate the main effects near the surface at large distances from the moving load 
when V > Cr. (Similar terms could be obtained when V < Cr but these terms are not 
needed in the present investigation.) The rationale for this approximate and simplified 
procedure is discussed in reference 12 for a time-dependent half-space problem and will 
not be discussed further. The results obtained in that paper are applied to a nuclear 
blast problem in reference 13. 


Equations for the Surface Displacements 

In order to obtain some indication of the magnitude and significant features of the 
displacement pattern in the three possible speed ranges, z was set equal to zero in the 
integrals for u", v", and w" and the results added to equation (50) to give the dis- 
placements at the surface of the elastic half -space. The equations for this special case 
are presented in this section. Some numerical results based on these equations are dis- 
cussed subsequently. 

The process of simplifying the integrals for u", v", and w" when z = 0 is 
straightforward. The main steps involved in evaluating the types of integrals encountered 
can be summarized as follows: 

(a) Integrals of rational functions of sin 0, cos 0, \Jm^- 2 sui 2 0 - 1, \Jm^ 2 sin20 - 1, 

\jl - M(- 2 sin20, and ^1- M^ ^sin^e can be expressed in terms of elliptic integrals and 

elementary functions by decomposing the integrand into partial fractions and then making 
use of the formulas and transformations given in reference 14. 
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(b) 


b b 

lim f f(6>)— — 5—— d0=7rC f(0)6(g(<?))d0 
z _0 + a z^+g^(0) J a 


(See ref. 11.) 
(c) 



f(0)6(g(0))d0 


can be expressed in terms of elementary functions by using equation (48). 

When the equations for the surface displacements are written, some simplification 
is attained by giving the radial and tangential components of the horizontal displacements 


ujr and u<p rather than the longitudinal and 
radial and tangential components are defined 
polar coordinates r, cp on the surface of the 
nates by the equations x = r cos cp and y = 



Figure 2.- Polar coordinate system f ,<p and polar co mponents up.u^ 
of horizontal surface displacements, r = 1 x 2 + y 2 ; tan <p = y/x. 


transverse components u and v. The 
as follows: Introduce a system of plane 
half -space related to the Cartesian coordi- 
r sin cp as shown in figure 2. The origin 
is thus always located at the instantane- 
ous position of the moving load. The 
angle cp is measured positive clockwise 
as seen from above the surface, cp = 0 
in the direction of motion, and r is the 
distance from the moving load. Then 

up- = u(x,y,0)cos cp + v(x,y,0)sin cp 


u cp = v(x,y,0)cos cp - u(x,y,0)sin cp 


The radial component ujr is positive when the displacement is outward along a line 
through the position of the load. The tangential component u^ is positive when the dis- 
placement is clockwise as seen from above the surface. 

The horizontal and vertical components of the surface displacements are: 

Subseismic case; 0 < V < ct, < M t < 1: 
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u 


n ( y 1 K 2 »«) + , ( 6 ~ 4 \fi) a l 2 M t 2 n / y 2 k 2 ,k) 

1 - sin 2 <p 

(6 + 4j3]o! 2 2 M t 2 . v 6K( ) (l - M t 2 sin 2 <p)(2 - M t 2 sin 2 <p) , 

+ ^ 5 “ n ( r 3 « 2 »«) - - + 3 A AA L u( V k 2 ,k 

1 - ao^Mt^sin ^cp \ ' sin^w sin" 4 a>T( cp) ' 


j ]M sin cp cos cp 

-|M t 2 

167T 2 M t 

L 

1 - “ Mt 2 sin2(£> 


o (2 - Mt 2 sin 2 <pj\ll - M, 2 sin 2 <puMt 2 sin 2 <p - 

H(-cos <p)sgn(sin cp )-* A» " 

167T T(<p) 


(55) 


w = 


P f J 3 


32nr |2T(<p) 
Mt cos cp 


(2 - M t 2 sin 2 <p) - M z 2 sin 2 (p + 8^1 - Mj 2 sin 2 <p)^l - M t 2 sin 2 <pH(-cos <p) 


^3 


1 + _JKUK l 


||l _ i m^ 2 1 " “ M|. 2 sin2(p - 1 1 - Q!j 2 M|. 2 sin 2 (p 


\l3]j3-5 


fl -a- 2 2 M t 2 1 - «2 2M t 2 sin 2 <p 


(56) 


Superseismic case; q. < c^ < V, 1 < < M t : 


u— — 0 
r 


3P* 

^ 167rr 


= - 


— H(-cos <p)sgn(sin cp) 


yl - Mj 2 sin 2 < pj | Mt 2 sin 2 <p - 1^2 - Mt 2 sin 2 <pj 
T (<p) 


(57) 


op* (2 - M|- sin 2 (p) l/l - M^"sin 2 ^c + 4(1 - M? sin 2 (pj\jl - M+ sin 2 <p 

w = — H(-cos <p)A ' 4 ^A Li (58) 

32 irr T(<p) 
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The symbols introduced in these equations which have not been previously 
employed are: 


k 2 2M t 2 
" 3 - M t 2 


K = 


1 

k 


V = 


2M^sin2^? 

3 - M^.2sin2<p 


T (cp) = ^1 - -- M^sin^^ ^1 - aj^Mj-^sin^^ ^ _ c^M^sin^) 


Y\ = " 8 

y 2 = 4(3|/3-5) 
y 3 = -4^3\|3 + 5^ 


K(k) is the complete elliptic integral of the first kind, n(o!,k) is the complete elliptic 
integral of the third kind, H(x) is the unit step function, and sgn(x) is the sign function. 

In equations (51) and (52), the terms having \jl - Mr2 in the denominator are 
included only when Mr < 1. The term in equation (53) which has - 1 in the 

denominator is included only when Mr > 1. The radicals ^M^sin^ - 1, 

^ 1 - M^sin and \fl - M^ ^sin2<p which occur in equations (55) to (58) are to be set 

equal to zero when the arguments become negative. 

It is interesting to note that except for several minor changes in notation, the equa- 
tions given for u- are identical in form to those obtained in reference 15 in connection * 
with a completely different half-space problem, the transient radial displacements pro- 
duced at the surface of an elastic half -space by a concentrated stationary load having the 
form of a step function in time. Moreover, when cp = 0, the expressions given for the 
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vertical displacements are also identical in form to the vertical displacement equations 
given in reference 15. 


Displacements Beneath the Line of Motion of the Load 
at Superseismic Speeds 

The equations for the displacements produced at points situated directly beneath 
the line of motion of the load, that is, in the plane y = 0, can, of course, also be obtained 
as special cases of the general equations (28) to (49). The transverse displacement v 
is zero. At subseismic and transeismic speeds, the inplane displacements u and w 
can be expressed as complicated combinations of elliptic integrals and elementary func- 
tions. At superseismic speeds, however, the inplane displacements can be expressed 
completely in terms of elementary functions. The equations for this latter case are pre- 
sented in this section. These results give some indication of the variation of the dis- 
placement pattern with depth at high speeds. 

As before, it is convenient to introduce a system of polar coordinates p,xp in the 
y = 0 plane related to the Cartesian coordinates by x = p cos xp and z = p sin xp, as 
shown in figure 3. The origin is again always located at the instantaneous position of the 

moving load. The angle xp is taken as 

zero along the positive X-axis and + ^ 

along the positive Z-axis. The radial and 
tangential components of displacement are 
found to be: 

u- = u(x,0,z)cos xp + w(x,0,z)sin xp 

> 

Up = ~ 1 ^ H (^ - §)sin xp cos xp 

(59) 



of the displacements in the vertical plane y = 0. p = t/x? + z 2 ; 
tan i/i = z/x. 


and 


u xp = w(x,0,z)cos xp - u(x,0,z)sin xp 

EL — 3 c °s .£ — H (V - -^[2 cos xp( 1 - M t 2 sin 2 xp\T$i(xp) +N2(i//)l[ 

p 32irM^sin^xp V 2 'L J| 


= - ; 


(60) 
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where 


NiM 




4n t ^cos ;//^cos2^ - 2M^ 2 sin2^ + ^2 cos 2 ^/ - M^sin^j ycos 2 ;// - 2M^^sin 2 ^/ 
~ j M^sin^^ ^nj-2 - aq^Mj-^shr^j _ a2^M^sm^i^j\jl - Mf-^sinfy 


n 2 G0 = 


4n^cos jI/^2 cos 2^ + ^sin.2^ ycos^i// + 2M^^sin.^ + ^2 cos2j// + M^sin^^ 

/n^2 _ i M^.2sin2i//'j ^2 _ aj2 M^sm^) ^ 2 _ o^M^sm^^ \jl - M^2 sin2^ 


( 61 ) 

3 


and where 


(62) 


= M^sin^i// - cos^i}/ 
n^2 = ^sin?i}/ - cos 2^/ 

The expression N i(if/) is to be set equal to zero when the argument of \Jl - M^sin^i// 

is negative; similarly, for N2 (i^) if the argument of \Jl - M^sb^i// is negative. The 

radial component u— is positive when the displacement is outward from the load. The 
tangential component u^ is positive when the displacement is clockwise as seen from 
the tip of the positive Y-axis. 

Numerical Results and Discussion 

Surface displacements .- Figures 4 to 10 show plots of the three components of the 
steady- state surface displacement for various values of Mt obtained from a numerical 
evaluation of equations (51) to (58). Recall that since Poisson's ratio is 1/4, the values 
Mt = 0.92, 1, and ]/3 correspond to load velocities equal to the speeds of the Rayleigh 
waves, shear waves, and compression waves, respectively. Hence, subseismic load 
speeds correspond to values of Mt less than unity, transeismic speeds to values of Mt 
between unity and ][3, and superseismic speeds to Mt greater than \Js. 

Since all components are directly proportional to P’ = P/ /z* and inversely pro- 
portional to the distance from the load r, it is convenient to introduce the nondimensional 

p? pt pr 

quantities ujr, w, and u^ defined by uj- = — u^, w = y=-w, and u^yu^. The 

quantities w and u^ are functions of Mt and <p only. Hence, for a fixed value of 
28 


Mt, they show the manner in which the vertical and tangential components w and u^ 
vary around the load. Values of w and u^ are shown only on the right-hand side of 
the line of motion of the load as seen looking in the direction of motion. The behavior of 
these components on the left-hand side of the line of motion can be obtained from the fig- 
ures by noting from equations (51) to (58). that, with respect to this line, the vertical 
component is symmetric w(-cc) = w (cp), and the tangential component is antisymmetric 

- u <p(-<p) = 


Figure 4 shows the radial component of displacement over the entire speed range. 
Only one plot is necessary since, surprisingly enough, uf is independent of position 

around the load. This 
result is contrary to intui- 
tion which would suggest 
that the radial displacement 
ahead of a moving load 
might be larger than that 
behind the load. The radial 
component is inward at sub- 
seismic speeds and becomes 
infinite as V increases to 
the R-wave speed. At all 
other speeds ujr is finite 
and continuous. The dis- 
placement is outward at the 
higher transeismic speeds 
u F = o for c r ^ v. and decreases to zero at the 

P-wave speed. The radial component of displacement is zero at super seismic speeds. 
This result too is rather unexpected. The value of uj- at Mt = 0 corresponds to the 
Boussinesq solution in the special case X' = p'. 



Figured- Radial displacement up = up at the surface. Entire speed range. 


' Figures 5 and 6 show the vertical and tangential components of displacement in the 
subseismic speed range. Figures 5(a) and 6(a) apply when V is less than the R-wave 
speed; figures 5(b) and 6(b) apply when V is between the R-wave and S-wave speeds. 

In addition to left-right symmetry, figure 5(a) shows that the vertical component is also 
symmetric fore and aft, that is, with respect to the line cp = 90°, provided V is less 
than the R-wave speed. Figures 6(a) and 6(b) show that the tangential component is anti- 
symmetric with respect to this line and zero along it at all subseismic speeds. The 
vertical component is seen to be considerably larger than either the radial or tangential 
components. The values of w and u ^ at Mt = 0 correspond to the Boussinesq 
solution for a static concentrated load when X’ = pt. 
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As can be seen from figures 5(a) and 5(b), the shape of the distribution of the verti- 
cal component around the load changes abruptly as V increases through the R-wave 
speed. This change is due to the fact that when V > c R , the load "outruns" the Rayleigh 
waves. A similar phenomenon, accompanied by certain changes in the displacement pat- 
tern, occurs when V increases through the S- and P-wave speeds. Some brief com- 
ments on this feature of the moving load problem seem to be in order before proceeding 
with a discussion of the figures. 

There are three types of waves which can propagate over the surface of the solid: 
The R (Rayleigh) waves, S (shear) waves, and P (compression) waves. Whenever the 
load speed exceeds the speed of one of these waves, the effect of the corresponding waves 
is confined to the region of the solid bounded by the surface and the lower half of a 
trailing Mach cone whose apex is at the load position. The axis of the cone is the line of 
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Figure 5.- Concluded. 


motion and the semiapex angle is sin"l( ^ aVe s P eec M w hich is equal to ft. or 0, 

\Load speed ) 4 l 

when the wave speed is equal to c^. or c^, respectively. The only region on the sur- 
face which is influenced by the waves is the trailing Mach wedge bounded by the two 

. One Mach wedge, formed by the R-waves, 

occurs at subseismic speeds when v>c R ; two wedges, formed by the R- and S-waves, 
occur at transeismic speeds; and three wedges are formed at superseismic speeds. 

The values of cp which give the position of the pairs of trailing Mach lines for the val- 
ues of Mj- shown in the figures are listed for reference in table I. The symbols cp-^, 
cp and cp^ are used for the angular positions of the Mach lines associated with the 
R-waves, S-waves, and P-waves, respectively. As indicated in the sketch in table I, 

<Pl < <Pt < <?R- 


Mach lines cp = ± 


7 t - sin' 


ifWave speed' 
Load speed 
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TABLE I.- LOCATION OF MACH LINES ON THE SURFACE 
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p — 

Figure 6.- Tangential displacement u® = at the surface. Subseismic speeds. 

In figures 5 and 6, one sees that as V increases, w and increase and 
become infinite at the R-wave speed. The magnitude of xl (0 is very small when V is 
between the R-wave and S-wave speeds. The figures show that across the R-wave Mach 
lines, VLy remains finite and continuous, but w experiences an infinite discontinuity 
and reverses its direction from up in front of the Mach lines to down behind the lines. 
This singularity in w across the R-wave Mach lines occurs whenever V is greater 
than the R-wave speed and behaves like the reciprocal of the perpendicular distance 
from the Mach lines. In this connection, recall that the radial component ujr is finite 
everywhere when V > c-^ and, like u^, remains continuous across the R-wave Mach 
lines. 
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Counterclockwise , Clockwise 





Figures 7 and 8 show the vertical and tangential components in the transeismic 
speed range. The largest vertical displacements occur near the R-wave Mach lines. 

The direction of displacement changes from up in front of the Mach lines to down behind 
them. The tangential component is finite and continuous on the R-wave Mach lines but 
is no longer antisymmetric with respect to the line <p = 90°. Both w and u^ are 
continuous on the S-wave Mach lines but change very rapidly close to these lines. (See 
table I for the angular position.) 

Some results at superseismic speeds are shown in figures 9 and 10. Except near 
the R-wave Mach lines, where w again becomes infinite and reverses direction, the 
vertical and tangential components are of comparable orders of magnitude. There is no 
disturbance ahead of the P-wave Mach lines. Both w and u^ change rapidly 
immediately behind these lines. Figure 10 shows that the tangential component acts only 
in the region between the P-wave and S-wave Mach lines and is zero everywhere else. 
Thus, since u^ is also zero, there are no horizontal surface motions behind the S-wave 
Mach lines. 
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Displacements ben eat h the load at superseis mic speeds.- The radial and tangential 
components of displacement in the vertical plane y = 0 are shown for superseismic 
speeds in figure 11. These components are also directly proportional to P* and 
inversely proportional to the distance from the load p. For a fixed value of M^, the 

functions u— = u- and u^ = tL£ plotted in the figures depend only upon xp and 

hence show the variation of the displacements around the load. As figure 3 indicates, 

\p = 90° directly beneath the load and \p = 180° on the surface directly behind the load. 
The radial and tangential components are seen to be of a comparable order of magnitude; 
in fact, the tangential component is larger than the radial component near the surface. 
Note that when <p = ip = n, u— = u— and u^ = -w. The positions of the trailing Mach 
lines can be determined from table I as before by simply reading \p in place of <p. 

Note that there is only one Mach line in the medium associated with each wave speed. 
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Figure 9.- Vertical displacement w = -^ w at Figure 10.- Tangential displacement Up-jpjU^ 

ihe surface. Superseismic speeds. at the surface. Superseismic speeds. 


Figure 11(a) shows that u- is zero ahead of the S-wave Mach line and is infinite 
immediately behind this line. This singularity behaves like the reciprocal square root 
of the perpendicular distance from the Mach line. The radial component is continuous 
across the R-wave Mach line and vanishes on the surface behind the load. This result 
agrees with the previous results which showed that there are no horizontal surface dis- 
placements between the S-wave Mach lines. Figure 11(b) shows that the tangential com- 
ponent is infinite behind the P-wave Mach line. This singularity is also of the reciprocal 
square-root type. At the S-wave Mach line, u^ is continuous but has a sharp discon- 
tinuity in slope; u^ is also continuous across the R-wave Mach lines. 

Since u— and u , are continuous across the R-wave Mach line, the vertical com- 
P 4* ’ 

ponent of displacement w = u^ cos 4/ + u- sin xjy is also continuous there. Consequently, 

at least at superseismic speeds, w is not singular everywhere on the R-wave Mach 

cone. Since the Rayleigh wave is a surface wave, it appears probable that the infinite 

discontinuity in the vertical component shown in the previous figures occurs only at the 

surface and that w is continuous within the half-space. 
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Figure 11.- Displacements in the vertical plane y = 0. 
Superseismic speeds. 


TIME -DEPENDENT HORIZONTAL SURFACE DISPLACEMENTS 

At the time the moving load is first applied to the surface of the half- space, tran- 
sient displacements are excited which are not accounted for in an analysis limited to the 
steady state. The complete time -dependent solution for the surface displacements, 
including these transient effects, is discussed in reference 5. In that paper the vertical 
displacements are obtained in closed form in terms of elementary functions over the 
entire speed range. The horizontal displacements, however, are very complicated. 
Expressions for these quantities are given in a form which is suitable for numerical cal- 
culation, but which bear no similarity to the combinations of elliptic integrals obtained 
in the present paper in the steady-state case. Since the steady-state solution is a 
limiting case of the time-dependent solution, it is reasonable to conjecture that the latter 
might also be expressible in terms of elliptic integrals. This conjecture was in fact 
found to be true. 

The purpose of this section is to obtain, from physical arguments, the basic equa- 
tions for determining the time-dependent horizontal components of the surface displace- 
ments and then to present for the various speed ranges expressions for these components 
which are comparable to the results obtained in the steady-state case. The final equa- 
tions presented are valid under the assumption used previously that Poisson's ratio is 
1/4. For the vertical displacements the reader is referred to reference 5. 

Superposition of the Response to an Impulse 

Introduce a rectangular Cartesian coordinate system X, Y, and Z fixed in the 
medium with the positive X-axis coincident with the line and direction of motion of the 
load and with the positive Z-axis pointing into the medium as shown in figure 12. This 
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system of axes is the same as that used 
in the investigation of the steady-state 
case. Now, however, it is further 
assumed that the origin O of the system 
is located at the point on the surface at 
which the moving load is initially applied. 
Further, define x = x - Vt, r ^ + y^, 

cp = tan"l ^L, r^ = + y2, and 

x 

d = tan _ l The physical meaning of 

these variables is indicated in figure 12. The time t = 0- corresponds to the time at 
which the load is first applied to the surface. The polar coordinate system r ,cp has 
its origin at the instantaneous location of the moving load as before; the origin of the 
r,0 system is the point at which the load is initially applied. The radial and tangential 
components of the horizontal displacements in this latter coordinate system are denoted 
by u r and ug, respectively. 

The manner in which one can obtain the surface displacements produced by a 
moving load and the way in which the solution of this problem can be made to depend upon 
the results obtained in reference 15 is as follows: The effect of a concentrated moving 
load may be approximated by placing a sequence of discrete impulses along the line 
traveled by the load and choosing the time and location at which each impulse is to act 

according to the desired velocity V of 
the moving load. Thus, as indicated in 
figure 13, if the load is initially applied 
to the surface at point O at time t = 0 
and the impulses are spaced at equal dis- 
tance Ax apart, the first impulse is 
applied at time t = 0; the second, at time 

t = the third, at time t = and so 

forth. The response at any point of the 
surface at any given time is then the sum of the responses to all the impulses which have 
been applied previous to that time. This approximation can be improved by increasing 
the number of discrete impulses until, in the limit of infinitely many impulses, the sum of 
the responses to the individual impulses passes over into a Duhamel superposition inte- 
gral representing the response to a continuously moving concentrated load. 



Figure 13.- Approximation of a moving load by discrete impulses. Time 
interval between impulses is At = 


Initial position Instantaneous 
of load position of load 



Figure 12.- Coordinate system and displacement components on the 
surface. 
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Position of 
impulse 



impulse. Au r = Rgcosfa - e); Au e = R5 sin(a - e). 


Let Rg(x - x',y - y’,t - t') be the 
radial displacement produced at point x,y 
at time t by a unit impulse acting normal 
to the surface at x' ,y' at time t' . (By 
symmetry, the tangential component of dis- 
placement is zero.) Then, as shown in fig- 
ure 14, the radial and tangential components 
Au r ,Au 0 of the horizontal displacement at 
point x,y due to a single unit impulse at 
x = Vt’, y = 0 are 


Au r = R 6 (x - Vt’,y,t - V)^ Vt ’) cos 6 - ± l s in 6 

\Jfr - Vt ') 2 + y 2 

> 

Au 0 = R 6 (x - Vt',y,t - P)SL- cos e-(x- Vt')sin 6 

y(x - Vt r ) 2 + y 2 


(63) 


If the moving concentrated load has a constant magnitude P, the strength of the approx- 
imating discrete impulses is P At. Hence, by superposition, one finds that the radial 
and tangential components of the horizontal displacement produced by the moving load 
are 


u r = 


U 0 = 



Au r dt* 

> 

Au q dt' 


(64) 


These equations are combinations of the equations given in reference 5 for the Cartesian 
components of displacement u x and u y . The equations given in the reference were 
obtained by means of the dynamic Betti- Rayleigh reciprocal theorem rather than from 
physical considerations. 

What is required then is an expression for the radial displacement produced by an 
impulse. This expression can be obtained from reference 15 which gives the expressions 
for the vertical and radial surface displacements produced by a stationary concentrated 
normal load applied to the surface as a step function in time. Since an impulse is the 
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time rate of change of a step function and the governing equations are linear, it follows 
that the response to an impulse can be, obtained by differentiating the results of refer- 
ence 15 for a step with respect to time. The analytical expressions for R5 are rather 
involved combinations of complete elliptic integrals of the first, second, and third kinds 
and are not given here. 

Figure 15 shows a plot of R§ for Poisson's ratio of 1/4 as a function of the 
c.t 

parameter r = in which c^ is the shear wave speed, t is the time after the 
application of the impulse, and r is the distance of the observer from the position of the 


impulse. The impulse is applied in the direction of the positive Z-axis so that the sur- 
rounding surface is instantaneously in compression. The response falls off as the 
reciprocal of the square of the distance r. At a fixed distance, R(t) is the time his- 
tory of the radial displacement. The symbols P, S, and R denote, respectively, the 


times of arrival of the compression wave 



shear wave (t = 1), and Rayleigh 


wave (t = oq ~ 1.09) generated by the impulse. As shown in figure 15, there is no 
response prior to the time of arrival of the P-wave. The displacement is outward 
immediately after the passage of the P-wave; then it diminishes rapidly and becomes 
inward until about the time of arrival of the S-wave. It then increases again and 
becomes logarithmically infinite at r = 1. After the passage of the shear wave, the 
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displacement rapidly decreases but remains outward. The displacement is infinite 
immediately after the passage of the Rayleigh wave and then steadily decreases to zero 
with increasing time. 

Expressions for the Time -Dependent Horizontal Components 
of Displacement at the Surface 

The results of carrying out the integrations in equations (65) are given in this sec- 
tion. The details are omitted; some pertinent remarks are given in the appendix in ref- 
erence 5. The expressions for the time-dependent horizontal components of displace- 
ment at the surface are: 

For the subseismic case; 0 < V < ct, < Mt < 1: 

Radial component of displacement: 


u r = 0 


u r = 
r /x'r 


P |f t(xx + y 2) 


16tt 2 


_ 2 i 18 2 2 

r z - j M t 2 y 2 v 


i±M nly s K 2 ,K) - 3 ^ 2?2 ' M t 2y2 )( ?2 ~ M t 2 y 2 ) n 

r 2 - a 2 2 Mt 2 y 2 ' 


J- Tk(xx + y 2 ) 


T(x,y) 


tjk z ,k 


u r = T7T 


18 


6-4^3 


^ r 16n 2 

6+4 fi U, 1 

2 2 2 2 (^ 3 ’ 

r z - afWlfy 1 




r 2 - — Mt 2 y 2 \ / r“ - cnj ‘ 


3(2r 2 - M t 2 y 2 )(r 2 - M t 2 y 2 ) 
T(x,y) 


n(77,k) 


T < 




9 6 ' 4 ^ n Uk^k) 

r 2 - T 2 ) 




1 < t < a. 


(65) 

Equation continued on next page 
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u r = 


U 0 = 


u 0 = 


u 0 = 


u 0 = 


P 

ju, f r 


H 


xx + y 2 j 


16 H 2 

6 +4\j3 

r 2 - o^M^y 2 
Pt 


18 


r 2 - i M t 2 y 2 


3 ^2r 2 - M t 2 y 2 ) ( r 2 - M t 2 y 2 ) 

n(yo,k) - L 

v 3 / T(x,y) 


r 2 - o!j 2 M^ 2 y 2 

11(77, k) 


> (65) 


_ o 
xx + y^ 1 


Qirii'r ^.2 _ a ^2^2y2^ T 2 _ o^2 


(“1 < T ) 


Tangential component of displacement: 



P \l Mty 


4 ' 16 it 2 


t 2 _ I 


18 _J n (y^ 2 ,,*) - (e - 4$ t2 " ai -n(y 2 K 2 ,«) 

?2.1 m 2 v 2 v ' F 2 - «! 2 M t 2 y2 ^ 


r“ - j M t y 


, _ x t 2 - a, 2 , n \ Ktr/tA (x 5 + y2) 2 3(2r 2 - M t 2 y 2 )(r 2 - M t 2 y 2 ) ^ J 

- ( 6 + 4 \fs) -j g 2T 11 ( r 3 K2,K ) + 0 2 2 t(x y) 11 ' 

' y r 2 - a 2 2 M t 2 y 2 ' 7 M t 2 y 2 Mt 2 y 2 r 2 T(x,y) 


kM t y 


^ 1677 2 


(< 5 + 4 #) 


r 2 -I 


18 


4 / \ - Ofi / 

nfy-ijk^ - (6 - 4 \/ 3 ) — n (r 2 > 

F-.lMtV 1 j 1 r 2 - a 1 2 Mt 2 y 2 V2 


. «« , 3 (!g - M <y.Kf - M|2? A „<„ 

\ 2 / m* 2 v2 ia. 2„2_2 T(x,y) 


r 2 - “ 2 2 

r 2 - or 2 2 M t 2 y 2 ^ 3 ’ ^ M t 2 y 2 M t 2 y 2 r 2 


kM t y 


f 4 ' 16 ji 2 


r2.I 


r2 - a, 2 


18 


— — n (yi> k ) - ( 6 - 4 fi) 2 — rn n ( r2,k ) 

. 1 M 2 v 2 V 1 > r2-ot 1 2 Mt 2 y 2 V / 


r - ~ M t y 


r 2 _ rv « 2 


- (e + 4N/s)_ 2 _ ff2 2 M 2 2y2 n ( y 3’ k ) Mt 2 y 2 ■ M t 2 y2r 2 


6K(k) + (xx + y 2 ) 3 (2r 2 - M t 2 y 2 )(r 2 - M t 2 y 2 ) 


T(x,y) 


PM^y r 2 - a x 2 
877/1’ p 2 - o' 1 2 M^ 2 y 2 


< r < 


( 66 ) 


1 < T < GJ-i 


( a i < T ) 
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For the transeismic case; c^. < V < c^, Mj < 1 < M^: 


Results for this case are obtained from the subseismic case by adding the following 
terms to the expressions for the region 1 j\j3 < r < 1: 

To u r , add 


3PM t r 

167rju T 



+ 



(2r 2 - Mt 2 y 2 )^r2 - Mj 2 y 2 ||Mt 2 y 


2 V 2 . r2 


T(x,y) 


(67a) 


To ug, add 


. W& sgn(y) „r + y2 l (”* - M t 2 y2)^7^V | /M t 2y2 • r2 

16 irju'r L J T(x,y) 


(67b) 


For the super seismic case; cj- < c^ < V, 1 < Mj < Mj-: 

Results for this case are the same as those for the transeismic case when r > 

f 

In addition, when t < 

ft 


u r = 


3PMtr, 
167T/Lt T ' 


'|H(-3 


^2r 2 - Mt 2 y 2 )^r 2 _ M^ 2 y 2 ||Mt 2 y 2 - r 2 
T(x,y) 


(68a) 


u 0 = - 


3P(kc + v2) ( 2r2 ' M t 2 y 2 )]R - . f2 

3P V a + y > H(-x)sgn(y)— 1 , * 

16irM'r T(x,y) 


(68b) 


In these equations 



k 2 = 


2 


3t 2 - 1 


k= 1/k 
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71 = 


2Mj 2 y2 

r^ - 


2M|. 2 sin %<p 
3 - M^ 2 sin^<^ 


T(x,y) = 





- «2 2M t 



7 \ = - 8 

72 = 4 ( 3 ^ 3 - 5 ) 
y 3 =-4(3^ + 5) 

The other symbols in these equations have been defined previously. The square roots in 
equations (67) and (68) are to be set equal to zero when the argument becomes negative. 

The values r = — , 1, and ai correspond, respectively, to the arrival times of the 

|/3 

P-wave, S-wave, and R-wave generated at the point of application of the load. Note 
that when V = 0, equations (66) give u q = 0, and equations (65) reduce to the result 
given in reference 15. 

The results obtained previously for the steady-state surface displacements may be 
extracted from equations (66) to (69) by forming the expressions for u— and u^, 
replacing x everywhere by x + Vt, and then taking the limit as t approaches °°. 

Discussion of Results 

The general features of the horizontal surface motions can be discussed with the 
aid of figure 16 which shows for comparison a sketch referring to each of the possible 
speed ranges. In each sketch, one is looking down upon the surface. The cross denotes 
the point at which the load was initially applied; the black dot indicates the position of 
the load. The dashed line indicates the path along which the load is moving. The con- 
centric circles whose common center lies at the cross indicate the position on the sur- 
face of the P-, S-, and R-waves generated at the time of application of the load. The 
lines from the position of the load tangent to the circles are the trailing Mach lines. 

In each of the four cases shown, the position of the load relative to the P-, S-, and 
R-waves does not change with time; for example, if Cr < V < c^, the load always lies 
between the R- and S-waves generated at the starting position. If the moving load is 
thought of as a sequence of impulses, it follows that distinct P-, S-, and R-waves are 
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X Point at which load initially applied 
• I nstantaneous position of load 

Path of moving load 

Mach lines 



(a) Subseismic speeds. 




(b) Transeismic speeds. (c) Superseismic speeds. 

Figure 16.- Waves and Mach lines associated with the surface displacements. 


generated by each impulse. But in the process of superimposing the impulses, these 
distinct waves except those generated at the origin lose their individuality and cannot be 
separated from the overall pattern of displacements. 

From equations (65) to (68), it follows that the radial component is symmetric and 
the tangential component is antisymmetric with respect to the path of the load. There 
is no disturbance ahead of the P-wave except in the superseismic case. Whenever the 
speed of the load is greater than one of the wave speeds, these waves will be confined to 
a wedge-shaped region between the Mach lines which trail behind the load. Three such 
wedges are formed in the superseismic case. It is convenient to distinguish between 
transient effects which tend to decay in amplitude as the load moves and steady-state 
effects which tend to form a fixed pattern of displacements moving over the surface at 
the velocity of the load. Of course, to an observer situated at a given point on the sur- 
face, the problem is completely transient in nature. 

The transient effects consist of the contribution to the displacements of the P-, 

S-, and R -waves generated at the point of application of the load. It can be inferred from 
the equations that both the radial and tangential components of displacement increase 
after the passage of the P-wave and remain continuous across the S-wave. The tangential 
component is finite during the passage of the Rayleigh wave but increases rapidly behind 
the wave. The radial component is finite until the arrival of the Rayleigh wave. Behind 
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the wave it becomes infinite with the characteristic reciprocal square root singularity. 
After the passage of the R-wave, the disturbance dies out and eventually leaves the sur- 
face in its original undeformed state. The displacements are generally larger in magni- 
tude in the direction of motion of the load than in the opposite direction. The amplitude 
of the P- and S-waves decays as 1/r; therefore, the associated displacements are not 
appreciable at moderate distances from the point of application. Also the Rayleigh wave, 
which has infinite amplitude in theory, will be dissipated by the losses present in a real 

medium. But since the amplitude falls off as — , the associated displacements should, 

\Jr 

at a given distance, be more perceptible than those associated with the P- and S-waves. 
For the same reason the Rayleigh wave should be noticeable at greater distances than 
P- and S-waves. 

Superimposed upon these transient displacements there is a steady-state pattern of 
displacements which has a fixed distribution with respect to an observer moving with the 
load and which moves over the surface with the load velocity. This pattern consists of 

the displacements produced by the presence of the concentrated load and, whenever 
V > cr, a system of one or more trailing Mach lines. The amplitude of the displacements 
in this pattern ultimately depends upon the distance from the load and not upon the dis- 
tance from the initial point of application. Thus, points near the path of the load receive 
the full effect of the steady-state pattern even though they are at considerable distances 
from the origin. 

In the steady-state pattern the displacements are proportional to 1/r and become 
infinite at the position of the load. In the superseismic case there is no disturbance 
ahead of the trailing Mach lines formed by the compression waves. At slower speeds 
only some, or none, of the waves are confined behind the load; in these cases the surface 
is displaced prior to the arrival of the load. The horizontal displacements vary contin- 
uously through the compression and shear wave Mach lines but change very rapidly as 
the lines are crossed. There is no discontinuity in the horizontal displacements across 
the Rayleigh wave Mach lines. A fuller discussion of the steady-state displacements and 
some numerical calculations of their magnitudes have already been presented. 

Figure 17 shows some sketches of the radial displacement along the positive X-axis 
shortly after the load has begun to move. These diagrams are based upon some calcula- 
tions made for specific values of t, c^, and V but are believed to be representative 
of the type of results to be expected in general. The sketches are not to scale as they 
are presented only to illustrate the gross features of the displacement pattern, in partic- 
ular, the overall shape and the direction of particle motion. It is interesting to note that 
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except for a multiplying factor which takes into account the motion of the load, the equa- 
tions for u r when y = 0 are identical to those given in reference 15 for a completely 
different problem, the radial displacement due to a step function. 

The symbols P, S, R, and L in figure 17 denote, respectively, the position of 
the P-wave, S-wave, R-wave, and the load. As indicated, the largest displacements 
occur near the location of the load and immediately behind the R-wave. Note that the 
infinite displacements behind the R-wave and near the load change direction as V 
varies. There is no radial displacement between L and P in the superseismic case. 





Figure 17.- Shape of radial displacements along positive X-axis. 


CONCLUDING REMARKS 

In this paper equations have been derived for the steady-state displacements pro- 
duced within an elastic half-space by a concentrated normal load in uniform motion over 
the surface. The results are valid for all depths, load speeds, and elastic constants. 
They may be used to evaluate the subsurface stresses and strains and, in conjunction 
with known superposition principles, could be employed to calculate the response of the 
half -space to distributed surface loads. 

The equations for the displacements produced on the surface at all speeds and in 
the vertical plane containing the line of motion of the loads at superseismic speeds have 
been expressed in terms of elementary functions and complete elliptic integrals in the 
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special case of Poisson’s ratio equal to 1/4. Some numerical calculations based upon 
these equations have been presented and discussed for a wide range of load speeds. 

On the surface it was found that the radial component of displacement was inde- 
pendent of the position around the load. This component is finite and continuous except 
at the Rayleigh wave speed (where all surface displacements become infinite) and van- 
ishes at superseismic speeds. At subseismic speeds the tangential component is anti- 
symmetric fore and aft and left and right of the load (quadripole distribution). At higher 
load speeds this component is finite and continuous across all Mach lines but changes 
rapidly near the compression wave (P) and shear wave (S) Mach lines. At superseismic 
speeds the only horizontal surface motions are tangential and these motions occur only 
between the P-wave and S-wave Mach lines. The vertical component predominates over 
the horizontal at all speeds. The most significant feature of the vertical component is 
the appearance of an infinite discontinuity across the Rayleigh wave Mach lines. 

At superseismic speeds the radial component of displacement in the vertical plane 
containing the line of motion of the load vanishes everywhere ahead of the shear wave 
Mach lines. This component is infinite immediately behind the S-wave Mach lines and 
vanishes with decreasing depth. The tangential component is infinite immediately behind 
the P-wave Mach lines and has a sharp discontinuity in slope at the S-wave Mach lines. 

In addition to these results for the steady-state problem, equations which include 
the transients excited by the initial application of the load have been derived for the 
horizontal surface displacements in the special case of Poisson’s ratio equal to 1/4. 
These equations are in a form which is directly comparable with that obtained for the 
steady-state condition and from which some of the steady-state results can be rederived 
by a straightforward limiting process. The new features in the displacement pattern 
which occur in the complete transient solution are the appearance of a distinct compres- 
sion, shear, and Rayleigh wave which propagate away from the point at which the load is 
initially applied. The disturbances carried by these waves decay according to the well- 
known laws governing spherical and cylindrical waves and the complete solution tends 
rather quickly toward the steady state. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 30, 1965. 



APPENDIX 


INVERSION OF THE FOURIER TRANSFORM OF THE DISPLACEMENTS 


The purpose of this appendix is to present in detail the steps leading from the gen- 
eral equations for the displacements (eqs. (24) to (26)) to the results quoted in equa- 
tions (27) to (46) and to the specific forms of the Rayleigh pole terms u', v’, and w' 
given for Poisson's ratio of 1/4 (eq. (49)). The investigation falls into two basic parts: 
(1) studying the singularities of the integrands for the purpose of determining how to 
define the values of the radicals and St and how to integrate around the Rayleigh 
poles to obtain a physically meaningful result, and (2) using this knowledge to carry out 
some of the integrations in closed form. 


Singularities in the Integrals for the Displacements 

The singularities in the integrals for the displacements consist of branch points 
associated with the P- and S-wave speeds and poles associated with the Rayleigh wave 
speed. When these singularities are discussed, it is convenient to allow A and p to 
take on complex values and to view the integrals in equations (24) to (26) as repeated 
contour integrals in two complex planes. The location of the poles and branch points 
relative to the path of integration then depends upon whether one considers the A-plane or 
the p -plane. In the following, p will be taken as a real number and the location of the 
singularities are treated with respect to the complex A-plane. This choice is arbitrary; 
one could equally as well take A real and study the location of the singularities in the 
complex p -plane. 

Rayleigh poles.- Any poles of the integrands occur at the zeros of A(A,p). When 
e = 0, the equation for these zeros is 


n2 


2 (\2 + p 2 ) _ M t 2 A 2 - 4 (a 2 + p 2 ) /a 2 + p 2 - M t 2 A 2 {a 2 + p 2 - Mj 2 A 2 = 0 
After this equation is rationalized and simplified, one obtains 

A 2 (a 2 + p 2 - ci' 1 2 M t 2 A 2 ) (a 2 + p 2 - | M t 2 A 2 VA 2 + p 2 - o^ 2 M t 2 A 2 ) = 0 


0 3 +\73 9 3 - J3 

where a-i* = — and a 0 ^ = — . The assumption that Poisson's ratio is 1/4 has 

1 4 ^ 4 

been used so that M^ 2 = ~ M^ 2 . The third and fourth factors may be disregarded as they 

O 

give extraneous roots. Since the numerators in the integrals for the displacements are 
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of the order of X 2 as X approaches 0, the double root X = 0 is associated with a 
removable singularity and may also be disregarded. 


The factor X 2 + ju 2 - aq 2 M t 2 X 2 = 0 results in two simple poles which are referred 
to as the Rayleigh poles since if one sets cr = c \/a\, Cr is the Rayleigh wave speed if 
Poisson's ratio is assumed to be 1/4. For the following discussion it is convenient to 


define 


M R=r-=«l M t and 


-R 



- 1 . 


When e > 0, the zeros of A(X,/i) are determined by the more complicated equa- 


tion [2 ^X 2 + m 2 ) - M t 2 X 2 - 2ieVxJ - 4^X 2 + ju 2 - ieVx) Sj St = 0. From this equation one 

can easily show by direct substitution that the Rayleigh poles are now given by 

\2 + ^2 = Mr 2 X 2 + ieVX. Hence, in the complex X-plane, the poles Xj and X 2 are 


near the points 



C'eV 


when Mr < 1 


and near the points 


0R 


C'eVi 


when Mr > 1. In these equations C' is a positive number. In the first instance 
Mr < 1, the poles are on the imaginary axis. However, when Mr > 1, the poles lie 
below the real X-axis for positive e and V and move onto the real axis from below 
as e approaches 0, as can be seen from figure 18. Thus, when e = 0 and Mr > 1, 
the path of integration along the real X-axis must be indented to pass above the two 
Rayleigh poles. (Note that if V is negative and, as a result, the direction of motion of 
the load is reversed, the Rayleigh poles lie above the real X-axis and the path of inte- 
gration must be indented below the poles. The proper method of indentation appears to 
be exceedingly difficult to discover without the use of the parameter e.) The contribu- 
tions of the indentations around the Rayleigh poles to the displacements are discussed 
subsequently. 


Imaginary 

X 


Path of integration,. 





Real 

X 


Figure 18.- Location of Rayleigh poles, X^ and X 2 , in the complex X-plane 
when Mr > 1. 
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Branch points ,,- Branch points occur where Sj = 0 and = 0. Their location is 
shown to depend upon the magnitude of V relative to c t and c^ . There are three 
possibilities: V < c^. < c^ , Cj. < V < c^ , and Cf- < c^ < V which are referred to as the 
subseismic case, transeismic case, and superseismic case, respectively. Since the 
argument of each of the radicals Sj- and S j is of the form A 2 + ju 2 - M 2 A 2 - iVeX 
(where M = in St and M = Mj in S^ , it is convenient to discuss features 
common to both radicals simultaneously. 


The branch points occur where \|x 2 + ju 2 - M 2 A 2 - iVeX = 0. Thus, in the complex 

T T 

X-plane, the branch points Aj and X 2 are near the points i 


M 


(_ )jl - M 2 


+ c'Ve 


when 


M < 1 but are near the points ± 


is a positive number. 




M 2 - 1 


■ - c f Vei when M > 1. In these equations c T 


Branch cut- 


Path of integration 


Imaginary 

X 

T 


Branch cut- 


i I M I 


vT-m 2 


Real 
* X 


1 \jl - Nr 


Figure 19.- Location of branch points Xj and X 2 in complex 
X-plane when M < L 


Imaginary 

X 


Path of integration 


.. -Ini 

V Vm 2 -1 





Real 


+i iw/! \Im 2 -( M 2 -i)x 2 i-iVS^y-ti' 2 


-Branch cut- 


1 1 . 

Figure 20.- Location of branch points Xj and X 2 in complex 
X-plane when M > 1. 


In the first instance M < 1, the 
branch points are on the imaginary axis 
as shown in figure 19. The path of inte- 
gration, the real X-axis, lies between the 
branch points, and the branch cuts must 
be chosen so as to lie on opposite sides 
of the real X-axis. It is convenient in 
this case to cut the plane along the imag- 
inary axis from x’i to +°°i and from 
X 2 to -°°i. When e = 0, 

|/a 2 + (U 2 - M 2 A 2 is taken as positive and 
and real everywhere on the path of 
integration. 

However, when M > 1, both branch 
points lie below the real X-axis for posi- 
tive e and V and, as indicated in fig- 
ure 20, move onto the real axis from 
below as e approaches 0. Now the path 
of integration lies above both branch 
points and the branch cuts must both lie 
in the lower half of the X-plane. When 
e = 0, it is convenient to cut the complex 
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A-plane parallel to the negative imaginary axis from x'i to Xj - °°i and from X 2 to 
X 2 - °°i. The branch of the radical ]jx 2 + ju 2 - M 2 X 2 is chosen which makes the radical 

positive and real when | x|< — ==L=r. The values of the radical on the path of integration 

]jM 2 - 1 

between and on either side of the branch points are shown in figure 20. (Notice that if 
the direction of motion of the load is reversed so that V is negative, the branch points 
and branch cuts lie above the path of integration. The correct procedure for positioning 
the branch cuts is not at all clear if the parameter e is omitted.) 

Since the treatment of the Rayleigh poles and the branch points has been decided, 
the parameter e has served its purpose and may be dropped. In the remaining work 
e=0. The basic steps in reducing the double infinite integrals for the displacements 
(eqs. (24), (25), and (26)) to single finite integrals plus Rayleigh pole terms will now be 
discussed. 


Reduction of the Integrals for the Displacements 

The reduction of the contour integrals for the displacements is carried out in two 
steps. First, the contribution from the indentation of the path above the Rayleigh poles 
when Mr > 1 will be extracted from the integrals. This step still leaves integrals 
of the same form as those in equations (24), (25), and (26) but which must now be under- 
stood in the sense of a Cauchy principal value. By transforming these integrals from the 
"rectangular" coordinates X,ju. to a system of "polar" coordinates r,0, it is found that 
the r-integration can be carried out explicitly and the remaining 9 integrals have 
finite limits. 

Contribution from the Rayleigh poles when Mr > 1.- Consider the integral 

I = C dx in which Xn is real and the contour C goes from -°° to +°° along 

d C * - *0 

the real X-axis but avoids the pole at Xq by an indentation above the real axis. Then 


I = lim 
17— 01 


£ 


Ao -T 7 f(X) 

-00 X - Xq 


dX + 


£ 


f(A) 
C 0 A “ *0 


dx+ £ 


f(A) 

X q+?7 A - Aq 


dX 


where Co is a semicircle of radius rj above the real axis, as shown in figure 21. On 
the semicircle, set X = Xq + Then 


lim C — dx = lim f ffxo + id# = -ni f(Xn) 

77 -O J C 0 A - X 0 77-0 J n V ' ‘ 
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Imaginary 

X 


p _ at t of _L , lt e 9_ r _ a ! ion 


c. 



— .Real 
X 


Figure 21- Indented contour. 


and 


lim 
77— 0 




x o-v f(x) r°° 


— dA + \ aAj = x 

A " A 0 ^ A o+?7 A " A 0 J ^-°° A " A 0 

where the integral on the right is a Cauchy principal value. Hence 


m 


•dA 


-l 


+°° f(X) 


dX 


r -m_ dx= .„ if (x 0 ) + r”jw_ (ix 

J C X - Xq ' ' 0_oo X - Xq 


(Al) 


By using equation (Al), the contribution from the indentation around the Rayleigh poles 
can be obtained as follows. 


P £C(X 11 ) 

Each of the A-integrals in equations (24), (25), and (26) is of the form \ —r 1 — -dA 

J c A (X,p.) 

-tM and x 2 = diil 
provided Mr >1. In the neighborhood of the pole at X = Xj this integral can be written 
J J^__dA where f(x) = and similarly, for the pole at A = X2. Applying 


where A(X,/i) has simple poles on the real X-axis at Xj = and X 2 = 


as 


C X - Ai 


A(A,/z) 


(x - X t ) 

equation (Al) to each of the poles gives 


I, 


g (X,m) 
c a(x,m) 


dX = -t 



— 

gU,M) 

g (x,m) 


771 

A(A , /J.) 

A(X,m) 



C\l 

1 

X=X2 A " x l 

A=Ai 


£ 


+x + “i<i4 dx 

,00 A(X^) 


= -771 


g(>i2,M-) [ g(^-l>M) 


9A(x 2 „u) 9A(X l3 jlt) 


9A 


9X J 


+ °° g( x »M) 
A(A,p) 


+ ^ dX 
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By applying this result to equations (24), (25), and (26), one finds that the expressions for 
the displacements become 

(u,v,w) = (u',v',w') + (u",v M ,w") 


where u", v M , and w M are the same integrals as those which appear in equations (24), 
(25), and (26) but are now taken in the sense of Cauchy principal values, and the additional 
terms u', v f , and w’ arising from the Rayleigh poles are of the form 



(*2>M) 

_ + A 



) BA 



dp 


BX 


BX 


By substituting into this integral the appropriate values of g and carrying out the indi- 
cated integration, the expressions given in equation (49) of the text are obtained. 


Transformation of the integrals for u", v", and w" to polar c oordinates .- The 
location of the branch cuts to be used in defining the value of a radical of the form 


|a 2 + p2 _ 


has already been discussed. 


Imaginary 

A 

t 


® Position of Rayleigh poles 
• Branch points 
— Branch cut 


® 

S t =\|B t 2 A 2 +|i2 

Sj = \|bj 2 A 2 + M 2 . 


I 



I 


t? 



Real 

A 


(a) Subseismic case, V* c t . 

Figure 22.- Position of branch cuts and values of Sj and on 
the path of integration in the three speed ranges. 


If these considerations are applied to 
each of the radicals and St (with 
e=0), the system of branch cuts and the 
values of St and S i to be used in the 
complex X -plane for evaluating the inte- 
grals u T1 , v", and w Tt are obtained. 
The system appropriate for each of the 
three speed ranges is shown in figure 22. 
The values of the radicals on the various 
portions of the real X-axis between and 
beyond the branch points are also shown 
in the figure. The ® marks indicate 
the positions of the Rayleigh poles with 
respect to the branch points when 
Mr >1. To abbreviate figure 22 some- 
what, it has be en con venient t o define 
B t ={l - M t 2 , ft - |/Mt 2 - 1, 

B i = fl - M Z 2 , and fy = 0^2 - 1. 
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Imaginary 

A 

4 <g> Position of Rayleigh poles 



(b) Transeismic case, c^<V<c^. 

Figure 22.- Continued. 


Imaginary 

A <8> Position of Rayleigh poles 



(c) Superseismic case, c f < V. 

Figure 22.- Concluded. 


If an integral such as u", v", and w" 


+00 

of the form ^ f(A,jU,St,Sj)dA d ju, 

— OO 


which so far has been treated as a repeated contour integral to be evaluated by integrating 
along the real A- and /i- axes, is now viewed as a double surface integral taken over a 
real A,j a plane in which A and p. are interpreted as a system of rectangular 
Cartesian coordinates, the branch points on the real A-axis in the complex A-plane 
become branch "boundaries" in the Cartesian A,p plane and the definitions of Sj and 
St, which formerly applied to various segments of the real A-axis, now hold over certain 
areas of the A,p plane. The value of a radical changes abruptly as one crosses a 
branch boundary. 
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Thus, for example, the imaginary branch points X = ±i which occur in the 

Bt 

subseismic case have no counterpart in the real X,p. plane. There are no branch 
boundaries and and St are positive and real everywhere. However, the real 


branch point X = + which occurs in the transeismic case forms two half lines in 

Pt I | 

one-half of the X,/n plane. Similarly, the branch point X = - jpl forms two half lines 

in the other half of the plane. The radical St is real provided | X | < and pure 

Pt 


imaginary if [ X | > These two conditions define four regions of the X,/i plane 

P t 

which are separated by the two branch boundaries X = ± 

ft 


The result of reinterpreting figures 22(b) and 22(c)'in this fashion is shown in fig- 
ures 23(a) and 23(b). The subseismic case offers no difficulty and is not shown since 
there are no branch boundaries and S^ and St are positive and real everywhere in the 
X,p plane. The two real branch points in the transeismic case give two branch bound- 
aries in the X,ju plane as shown in figure 23(a). These lines form angles = sin"-*- — 

Mt 

with the /Lt-axis. The values of St in the four wedge-shaped regions bordered by the 
branch boundaries are also shown in the figure. In this case, Sj is positive and real 
everywhere. As can be seen from figure 23(b), the four real branch points in the super- 
seismic case give four branch boundaries in the X,ju plane. Two of these lines are the 
same as those that occur in the transeismic case. The other two make angles 

= sin-1 _T_ with the ju- axis. The values of and St in the eight wedge-shaped 
regions of the plane bordered by the branch boundaries are also shown. Although they 


-Branch boundaries 






,2 \v 








^ H ^ . 

2,2 2 


r * \K * 


111 i 

(a) Transeismic case; \ +\i everywhere. 

Figure 23.- Position of branch boundaries and values of Sj and S\ in the plane. 
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Branch boundaries 
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+ i^/p t 2 X 2 -p 2 
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(b) Superseismic case. 


Figure 23.- Concluded. 


are not indicated in the figure, the Rayleigh poles also show up as two lines (when 

Mr > 1) making angles 6 >r = sin~l < 9\_ with the p-axis. The integrand f^X,p,St,S£ 

is singular everywhere along these lines; this property is signified by the Cauchy prin- 
cipal value. Another description of the connection between figures 22 and 23 may be of 
some value to the reader: figures 22(b) and 22(c) are cross sections of figures 23(a) 
and 23(b) along the lines p = Constant. 


Since the values of the radicals S(; and have been determined throughout the 
X,p plane, the values of the integrands of the integrals for u", v", and w" are also 
known and some progress can be made toward carrying out part of the remaining integra 
tions in closed form. For this purpose it is convenient to transform the integrals from 
X ,p coordinates to plane polar coordinates r ,Q defined by X = r sin 9, p = r cos 9. 
By using figure 23 as a guide, one finds that an integral of the type 


+ CO 

J£f(x,p,St,Sj)dXdp 

-00 


can be cast into the following form: 
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For the subseismic case: 


y °° V /2 

X fi(r,0)d0 r dr 
0 J 0 


For the transeismic case: 


r °° n 0+ p ff/2 

* fi(r,0)d0 + X f 2 (r,0)de 

0 J 0 J 0 t 


r dr 


For the superseismic case: 

6t 


f\ °° r> n Oj p tt /2 

\ \ fi(r,0)d0 + \ f 2 (r,0)d0 + \ f 3 (r, 0 )d 0 

J 0 J 0 J 0 t J 0 Z 


r dr 


in which 


f 1 (r,0) = ffr sin 8, r cos 8, r \jl - Mt 2 sin 2 0 ? r \Jl - M^ 2 sin 2 0 ) 

' \ 

+ f^r sin 8,-r cos 8,r\Jl - Mt 2 sin 2 0,ryl 

- M^ 2 sin 2 0j 

+ f^-r sin 0,-r cos &,r\jl - Mf- 2 sin20,r ]j 

'l - Mj 2 sin 2 0) 

+ f^-r sin 8,r cos 0,r^l - Mt 2 sin 2 0,r^l 

- M^ 2 sin20^ 
\ 

f 2 (r,0) =f^r sin 0,r cos 0 ,-ir^Mt 2 sin 20 - l,r 

- 2 sin 2 0 j 

\ 

+ f^r sin 0,-r cos 0 ,-ir^Mt 2 sin 20 - l,r 

\J 1 - M^ 2 sin20j 

\ 

+ f ^-r sin 0,-r cos 0 ,+ir^Mt 2 sin 20 - 1, 

r\Jl - M z 2 sin 2 0j 
\ 

+ f^-r sin 0,+r cos 0 ,+iryMt 2 sin 20 - 1, 

r\jl - Mi 2 sin 2 0j 

\ 

f 3 (r,0) =f^r sin 0,r cos 0 ,-ir^Mt 2 sin 20 - 1,-ir \j 

2 sin 2 0 - 1 J 

\ 

+ f^r sin 0,-r cos 0,-ir^M^ 2 sin 2 0 - 1,- 

ir^Mj 2 sin 2 0 - 1 j 

\ 

+ ff-r sin 8,- r cos 0,+ir^Mt 2 sin 2 0 - 1,+ir^Mj 2 sin 2 0 - 1 j 

+ f ^-r sin 0,r cos 0 ,+ir\Mt 2 sin 20 - 1,+ 

•ir^M^ 2 sin 2 0 - lj 


(A2a) 


(A2b) 


(A2c) 


(A3) 
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Equations (A2) and (A3) are the basic formulas for reducing the integrals for u", v", 
and w", to single finite integrals. If these formulas are applied to equations (24), (25), 
and (26), it is found that the r-integration can be carried out explicitly and the results can 
be expressed in terms of the combinations of elementary functions defined in equation (47). 
The final results are given in equations (28) to (47) of the text. Note that in deriving 
these equations there is no need to assume any specific relation between the wave speeds 
c^ and Cf- or the Lame constants X' and p'. 
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